Measurement of the hydrodynamic forces between two polymer-coated spheres 
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The hydrodynamic forces between Brownian spheres are determined from a measurement of 
the correlated thermal fluctuations in particle position using a new method, two-particle cross- 
correlation spectroscopy (TCS). A pair of 1.3 /im diameter polymer-coated poly(methyl methacry- 
late) were held at separations of between 2/im and 20^tm using optical traps. The mobility tensor is 
determined directly from the statistically-averaged Brownian fluctuations of the two spheres. The 
observed distance dependence of the mobility tensor is in quantitative agreement with low-Reynolds 
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I. INTRODUCTION 



The hydrodynamic interactions between colloidal particles are important from both a fundamental and an industrial 
viewpoint. They determine, for instance, the rheological behaviour of suspensions, the kinetics of aggregation and 
phase separation and many other common colloidal phenomena Q. Yet despite this, the hydrodynamic properties of 
all but the simplest colloidal system have been a subject of considerable debate (2). A key factor in this uncertainty has 
been the intrinsically long-ranged nature of the hydrodynamic coupling between solid particles. In the dilute limit, 
where it is sufficient to calculate just the leading order interaction, the solution of the stationary Stokes equation 
H reveals that the interactions decay like the inverse separation (the Oseen tensor) . The non-local nature of such 
interactions has led to considerable theoretical and numerical difficulties. Experiments have, conversely, also been 
1 problematic. Methods such as light scattering which have been used extensively in the past to probe the dynamics of 
concentration fluctuations provide only limited, indirect information on the microscopic nature of hydrodynamics in 
q ■ suspensions. 

In this paper we report a detailed experimental study of the distance dependence of the hydrodynamic interactions 
between an individual pair of polymer-coated colloidal particles. Optical tweezers are used to hold two uncharged 
I ' poly(methyl methacrylate) (PMMA) spheres apart at separations from between 2fim and 20/im. The hydrodynamic 
, forces are measured using a new, highly sensitive, experimental technique, two-particle cross-correlation spectroscopy 
t-H ' (TCS). TCS measures the statistically-averaged instantaneous fluctuations in position of two probe spheres. The 
in-plane position of each sphere is measured to nanometer precision using a quadrant photodetector. This technique 
is applied to a pair of colloidal particles suspended in a Newtonian liquid to yield a detailed direct test of low- Reynolds 
^1 ' number predictions of hydrodynamic forces in a simple colloidal system. 

Several reports of the test of theoretical predictions for the hydrodynamic coupling between a pair of spheres have 
already been published ^, Our experiments differ in several regards. First, in the experiments reported to date, the 
surfaces of the colloidal spheres have been bare and not covered by a polymer layer. Since the adsorption or anchoring 
of a polymer onto the surface of a particle is a common method of imparting colloidal stability it is important to 
establish if the presence of a polymer layer modifies the hydrodynamic forces. Indeed theoretical calculations || 
predict that, flow within the polymer layer, removes the divergence of the hydrodynamic forces seen at small pair- 
separations. Second, the reported experiments have used charged spheres at relatively low electrolyte concentrations 
(e.g. O.lmM in Q) so that residual electrostatic interactions between the charged spheres or between the confined 
spheres and the walls complicates the interpretation. The experiments reported here used an uncharged colloid, for 
which previous work [Q, |^] has shown that the interaction potential is well approximated by that of hard spheres. 
Our paper is organised as follows: In the next section we describe the details of our experiment. Section III details 
k> , the Brownian motion of an isolated sphere in a harmonic optical potential and describes the motion of a pair of 



dynamically-coupled spheres. We present our results in §|fv| before concluding. 



II. EXPERIMENTAL METHODS 



Measurements were performed on a dilute suspension of uncharged poly(methyl methacrylate) (PMMA) spheres of 
radius 0.65 ± 0.02/im. To minimize the van der Waals forces between the PMMA spheres the surface of each sphere 
was covered with a covalently bound polymer brush (~ 100A thick) of poly(12-hydroxy stearic acid) A dilute 
suspension of the spheres (volume fraction <f> ~ 10~ 7 ) in a mixture of cyclohexane and cis-decalin was confined within 
a rectangular glass capillary, 170/im thick. The ends of the capillary were hermetically scaled with an epoxy resin to 
prevent evaporation and to minimize fluid flow. A small concentration of free polymer stabilizer was added to the 
suspension to reduce the adsorption of the spheres onto the glass surfaces of the cell. 
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A pair of spheres were trapped in a plane about 40/um above the lower glass surface of the cell using two optical 
traps. The traps were created by focusing orthogonally polarised beams from a Nd:YAG laser (7910-Y4-106, Spectra 
Physics) with a wavelength of A = 1064nm to diffraction-limited spots using an oil-immersion microscope objective 
(100x/1.3 NA Plan Neofluar, Zeiss). The resulting optical gradient forces localise the sphere near the focus of 
the beam. The centre-to-centre separation of the two traps could be varied continuously from between 2 and 30/im. 
Accurate positions for the optical traps were determined by digitising an image of two trapped spheres with a Imaging 
Technology MFG-3M-V frame grabber. The spheres' locations were measured to within 40nm using a centroid tracking 
algorithm. 

While the mean position of each sphere is fixed by the position of the corresponding laser beam, fluctuating 
thermal forces cause small but continuous displacements of the particle away from the centre of the trap. For small 
displacements the optical trapping potential is accurately described by a harmonic potential JlO| . The restoring force 
on the particle is proportional to the displacement with a force constant which for a given particle and beam profile 
is a linear function of the laser power. In the experiments detailed below the intensities of the orthogonal beams were 
carefully adjusted until the stiffness of the two traps differed by less than 5%. The trap stiffness k was typically of 
the order of 5.1 x 10 _6 Nm _1 which corresponds to a rms displacement within the trap of about 40mn. The intensity 
of each beam at the focal plane was estimated to be of the order of 30mW. 

The positions of the two trapped spheres, r\ and r 2 , were tracked with nanometer resolution by observing the 
interference between the transmitted and scattered light in the back-focal plane of the microscope condenser using a 
pair of quadrant detectors [j^TJ . Difference voltages from the sum of the horizontal (X) and vertical (Y) halves of the 
quadrant detectors are linearly proportional to the displacement of the sphere from the optical axis of the trap. The 
trajectories, ri(t) and r 2 (t), were measured for a pair of spheres with mean separations r = \ri — r 2 \ from between 
2.5^m and 20/im. For each value of r, the Brownian motion of the two spheres was followed for a total of 420s, at 
intervals of 50^sec, to yield 2 23 (8.4 x 10 6 ) samples of the spheres' dynamics. 



III. BROWNIAN MOTION OF CONFINED SPHERES 



A. An isolated sphere 

A single hard sphere of radius a, moving with a constant velocity U through an unbounded fluid of viscosity r\ 
experiences a hydrodynamic drag force Fp in the direction opposite to motion. In the low Reynolds number limit ||] 
the velocity of the sphere is a linear function of the force exerted on the particle by the fluid, 

U=-b F D . (1) 

The constant bo is the mobility of an free particle which, if there is no slip at the boundary of the particle, is given 
by Stokes Law as 

b » = ( 2 ) 

When the Brownian sphere is confined by a potential, U (x) , the drag force increases. The total force on the particle 
consists of a random Gaussian force f(t) together with an additional force due to the potential field, —dU/dx. For a 
harmonic potential (of stiffness k) in the 'long-time' limit, where inertial terms are negligible, the motion of a confined 
Brownian sphere is described by the Langevin equation, 

^= bo [ m -kx(t)] (3) 
with a random particle force f(t) which is Gaussian distributed with the moments, 

</(*)> = o 

(/(*)/(*')) = 2bo 1 k B T6(t-t'). (4) 

This Langevin equation is readily solved by standard methods [^| and the position autocorrelation function (x(t)x(0)) 
determined. Since there is only one characteristic timescale in the long-time limit, the autocorrelation decays expo- 
nentially 

( ar (t)x(0)) = ^exp(-t/r) (5) 
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with a decay time r which is physically just the time taken by a sphere to diffuse a distance I*, 

where Z* is the classical turning point of the confining potential, or the separation at which the potential energy of 
the trapped sphere equals its thermal energy ksT. 



B. A pair of spheres 



The motion of a pair of harmonically-bound particles differs from III A because the hydrodynamic forces couple 
the motion of the two spheres together. As one particle moves a flow is created in the surrounding fluid which drives 
fluctuations in the position of a neighbouring second sphere || . In this section we analyze the correlated motion of 
a pair of particles which are coupled by such dynamic forces. We assume, for simplicity, that (1) the two trapped 
particles have the same diameters and are contained within optical traps with identical force constants; and (2) that 
there is no potential coupling between the two spheres. 

The hydrodynamic forces acting between equal-sized spheres have been calculated by a number of authors jll], [u|, 
jL5j. In low- Reynolds number flow the hydrodynamic interactions between two spheres can be described by a set of 
linear relations between the force or torque exerted on a sphere and the corresponding translational and rotational 
velocities. If, as here, the spheres are freely rotating the applied torque must be zero and one can eliminate the 
angular velocities. In this case the linear relation between the forces and translational velocities defines the mobility 
tensor 6, 



UA (bu b 12 \ (Fx 
U 2 I 621 b 22 ) \F 2 



(7) 



where the two spheres are labeled 1 and 2 and the equivalence of the two particles implies that b\\ = b 22 and 612 =621- 
The spherical symmetry of the Stokes limit reduces to one of axial symmetry and the mobility tensor depends crucially 
on the geometry of the two spheres p3| . The linearity of the Stokes equation implies that each of the matrices b^ 
can be decomposed into a pair of mobility coefficients which describes motion either along the line of the centres or 
perpendicular to it, 

TV VV 

b ij {v)=A ij {r)^ + B ij {l--^) (8) 

where the coefficients, Ay and Bij, detail the longitudinal and transverse mobilities, respectively. In the remainder of 
this paper we shall confine our discussion to the longitudinal motion alone. In this case only two mobility coefficients 
are needed to quantify the hydrodynamic forces. Each of these mobilities, An and A12, is itself a function only of the 
normalised separation p = r/a of the two spheres. 

The mobility coefficient, A 12, couples the fluctuations along the line of the centres of the two spheres, chosen here 
as the x-axis, so that the particle coordinates, x\ and x 2 , are no longer independent. The extent of their correlation 
is determined from the solution of the Langevin equation 



xi\ (A u A 12 \ (h{t)-kxx 
x 2 ) U12 An) \f 2 (t)-kx 2 



(9) 



with the random force characterized by the moments 



(/*(*)) = 

= 2(A-%k B TS(t-t'). (10) 

and (A _1 )jj as the inverse matrix of Aij. Note that the mobility gradient terms, ksT /2^2^ dAij /dxj, in the conven- 
tional Langevin equation have been ignored in equation ^ because in our experiments the positional fluctuations are 
typically two orders of magnitude smaller than the mean particle separation. 

To solve this coupled Langevin equation we introduce the normal coordinates Xi, 
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and choose the coefficients Cy so that the equation of motion for Xi has the following form 16 

dXj 
dt 



-kXiXi + Fiit) 



(12) 



with i = (1 or 2). It is readily shown that the matrix Cy consists of the normalised eigenvectors of the mobility matrix 
Aaa so that the normal modes are 



X x = 
X 2 = 



1 
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(X\ + X 2 ) 



X 2 



(13) 



which describe, in turn, a symmetric collective motion (Xi) of the centre of mass of the two spheres and an anti- 
symmetric relative motion (X 2 ) of the two spheres with respect to each other along the line of their centres. The 
mobilities A^ of the two modes are the eigenvalues of the matrix cy , 

Ai = A n + A 12 

A 2 = A n -A 12 . (14) 

while the Fi's are random forces which satisfy 

(m)) = o 

(Fi(t)Fj(t')) = 2S lJ X l k B T5(t-t'). (15) 

Since the random forces are independent of each other, motion of the two normal modes are also independent of each 
other. The hydrodynamic term which couples the motion of the two spheres, A 12 , leads to an asymmetry in the 
decay times of the normal modes. The time correlation functions of the normal coordinates are calculated readily 
from equation [ij as 

k B T 



( Xl (t)X 3 (0))-^ — 



exp(-t/r l ) 



with decay times 



n = 



(16) 



(17) 



Inverting the coordinate transformation of equation ^ gives the normalised time correlation functions of the particle 
centres 



hij(t) 



(Xi(t)Xj(p)) 



(18) 



as 



h n (t) 
h 12 {t) 



exp(-t/ri) + exp(-t/r 2 )] 
exp(-t/n) - exp(-t/r 2 )] . 



(19) 



Inspection reveals that the cross-correlation is very sensitive to the hydrodynamic coupling between the two spheres. 
At small times, t — * 0, h\ 2 records only the time-averaged, static correlations |l7j which, at thermal equilibrium, 
depends only on the interparticle potential. In the current experiments there is no potential coupling between the two 
spheres and so hi 2 (0) = 0. At long delay times, t — > oo, the hydrodynamic flows which couples the motion decay to 
zero so that the positions of the two spheres are uncorrelated and hi 2 (t —> oo) = 0. The cross correlation is therefore 
zero both at short and long times. Equation [l9] reveals that the cross correlation will also be zero at intermediate 
times unless the two decay times, t% and t 2 , differ. From equations [l7] and [l4| this difference is a linear function of 
the hydrodynamic coupling term A\ 2 . In the limit where A\ 2 <C An, which is the case in most physical situations, 
the cross-correlation has a minimum at a time which is fixed, to leading order, by the diagonal mobility An 
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while the depth of the minimum is determined by the ratio of the off-diagonal and diagonal mobilities 
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IV. RESULTS 



The normalised longitudinal position cross correlation, h^lt), was measured for a pair of trapped spheres over a wide 
range of separations (2.5 fim < r < 20 /im). Figure |l| shows typical data for four sphere separations. Three features in 
the experimental data are striking. First, the data show, rather surprisingly, that hydrodynamic interactions cause the 
two particles to be anti- correlated at intermediate times; second, the time at which the two particles are most strongly 
anti-correlated, the time t* at the minimum of ft, 12 , does not vary with the sphere separation; while conversely, the 
strength of the anti-correlation increases markedly as the separation r reduces. 

To interpret these observations the decay times t\ and r 2 of the symmetric and anti-symmetric normal modes were 



extracted from a least squares fit of the data to equation 19. The quality of the resulting fit may be gauged, for the 
four separation depicted in figure [l], by studying the solid curves, which are seen to accurately reproduce the measured 
data. From the experimentally determined decay times and trap stiffness the elements of the mobility tensor may be 
estimated as, 
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A 12 = ir[ —-— ). (22) 



The experimentally-determined scaled mobility elements, b^ 1 Aij{r), are plotted in figure as a function of the 
dimensionless separation of the two spheres, p — r/a. The two mobilities show a strikingly different dependence on 
the sphere separation. While the diagonal mobility is largely unaffected by the sphere separation the off-diagonal 
term scales approximately inversely with p. These trends are, of course, consistent with the observations made above 
that the position, t = t* , of the minimum in the correlation function does not shift with separation (see equation |20| ) 
while the depth increases with reducing separation (equation . 

The experimental values for the mobilities may be compared with theoretical predictions for the hydrodynamic 
coupling of two hard spheres. Batchelor [ jl3| , for instance, has given the following expressions for the longitudinal 
mobilities, 

K l Au = i-^ + o(p- 6 ) 

&o%2 = f -± + 0(p- 7 ) (23) 
zp p 

which are exact in the limit of large centre-to-centre separation p, as has been confirmed by [fl4j| . The solid curves in 
figure ^ show the predictions of the Batchelor theory for the longitudinal coupling. As is clear from this figure, the 
measured mobilities agree very well with the theoretical predictions over the entire experimentally accessible range of 
separations. 

The deviations between theory and experiment evident in figure ^| are due largely to the experimental difficulty 
measuring the force constant k. This is seen in figure ^ where the experimentally determined mobility ratio, Ai 2 /An, 
which from equation ^2] does not require any knowledge of k, is plotted as a function of the inverse separation 1/ p. The 
almost quantitative agreement seen, with no adjustable parameters, between the data and theory confirms the accuracy 
of Batchelor's theoretical description of pair hydrodynamics. In addition the close agreement between experiment 
and theory suggests that, at least for the range of distances explored in the current experiments, the hydrodynamics 
forces between polymer-coated and uncoated spheres are very comparable. 



V. DISCUSSION 



We have presented a detailed experimental study of hydrodynamic coupling between an isolated pair of polymer- 
coated hard-sphere colloids. We find near quantitative agreement with Low Reynolds-number predictions for the 
hydrodynamic coupling between a pair of spheres. Surprisingly, we observe a strong anti- correlation in the positions 
of the two coupled spheres at intermediate times. At first sight this result looks counter-intuitive since one might 
expect naively a symmetric correlation between spheres. However the effect is a dynamic time-dependent phenomenon. 
The origin of which may be understood from the normal modes of the system. The motion of two spheres, along 
the line of their centres, decouples when analyzed in terms of a symmetric collective mode and an anti-symmetric 
relative mode. The mobilities of these independent modes are, from the asymptotic expressions of Batchelor |13j and 



6 



equation [l|, 

Az = 6 »{ 1 -^ + 7-^ + ("" a '}' < 24 > 

where p is the dimensionless centre-to-centre separation, p = r/a. Examination of the leading terms in this equation 
reveals that the mobility Ai of the symmetric mode is enhanced and the anti-symmetric mode reduced when compared 
with an isolated particle. The reduction in mobility of the anti-symmetric mode reflects the difficulty of squeezing 
fluid out of or into the narrow gap between two approaching spheres while the increased mobility for the symmetric 
mode is caused by the tendency for the fluid flow generated by one sphere to entrain a neighbouring sphere. 

The asymmetry in the mobilities of the normal modes, seen in equation ^4|, causes the decay times for thermal 
fluctuation in the two modes to differ. When the two spheres are close together, the mobility of the symmetric mode 
is enhanced compared with the anti-symmetric mode. As a consequence, symmetric fluctuations decay more rapidly 
than their anti-symmetric counterparts. At t = the proportions of thermally-excited symmetric and anti-symmetric 
fluctuations are equal since the positions of the two spheres are uncorrelated. With increasing time the amplitudes 
of both fluctuations decay. However the anti-symmetric fluctuations decay at a slower rate than the symmetric 
fluctuations so that the cross correlation develops a pronounced anti-correlation. The anti-correlation is however 
dynamical since over long time all fluctuations decay and the spheres become again uncorrelated. 

In summary, we have shown that 2-particle cross-correlation spectroscopy (TCS) is a promising new technique for 
the quantitative determination of hydrodynamic interactions. TCS experiments are very flexible; the particle size, 
separation, potential interactions and indeed the dispersion medium can all be changed independently of each other. 
Variations of the methods described in this paper could, for instance, be used to follow the time course of collective 
fluctuations in a dense (host) complex fluid from the real-space trajectories of inserted probe colloidal particles. 
Currently, we are utilising TCS to study the many-body hydrodynamic interactions in concentrated particulate 
suspensions. Two probe PMMA particle are trapped within an index-matched silica suspension of volume fraction 
4>. The resulting fluctuations in the trajectories of the two probe PMMA particles are used to determine the effective 
pair-mobility tensor in the suspension, as a function of particle separation and <f>. These measurements promise 
to provide new and detailed experimental information on the spatial and temporal development of hydrodynamic 
interactions in concentrated suspensions. 
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FIG. 1: Cross correlation functions for 1.30 /J,m diameter PMMA spheres as a function of delay time t and at four centre-to- 
centre separation r. (a) r — 2.47/im, (b) r = 4.19/im, (c) r = 7.44/im, and (d) r = 11.47pm. The motion is measured parallel 
to the separation vector. The solid line shows a fit to equation [13. For clarity only every twelve data point is plotted. 
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FIG. 2: The experimentally-determined mobility coefficients, 6q 1 A i j, for motion along the line of centres as a function of the 
inverse centre-to-centre separation p = r/a . Here bo is the mobility of the PMMA particle of radius a. Solid lines show 
the predictions of low- Reynolds number hydrodynamic calculations for the case of two interacting solid spheres with no 
adjustable parameters. 
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